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Abstract 

We study the initial-boundary value problem for an Euler-Bernoulli beam model with 
discontinuous bending stifltness laying on a viscoelastic foundation and subjected to an axial 
force and an external load both of Dirac-type. The corresponding model equation is fourth 
order partial differential equation and involves discontinuous and distributional coefficients 
as well as a distributional right-hand side. Moreover the viscoelastic foundation is of Zener 
type and described by a fractional differential equation with respect to time. We show how 
functional analytic methods for abstract variational problems can be applied in combination 
with regularization techniques to prove existence and uniqueness of generalized solutions. 
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1 Introduction and preliminaries 

We study existence and uniqueness of a generalized solution to the initial-boundary value problem 

dfu + Q{t,x,d:,)u + g = h, (1) 
D^u + u = eD^g + g, (2) 
w|t=o = /i, dtu\t=Q = f2, (IC) 

u\x=0 = U\a:=l = 0, 9j;U|a;=0 = dxU\x = l = 0, (BC) 

where Q is a differential operator of the form 

Qu :— dl{c{x)dlu) + b{x, t)dlu, 

b,c,g,h, fi and /2 are generalized functions, 6 a constant, < 6* < 1, and D" denotes the left 
Riemann-Liouville fractional derivative of order a with respect to t. Problem (IT])-(l2]) is equivalent 
to 

d'^u + Q{t,x,dx)u + Lu = h, (3) 
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with L being the (convolution) operator given by (£ denoting the Laplace transform) 

Lu{x,t) = C-^(^^±^^{t)Hu(x,t), (4) 

with the same initial (IC) and boundary (BC) conditions (cf. Section [3]). 

The precise structure of the above problem is motivated by a model from mechanics describ- 
ing the displacement of a beam under axial and transversal forces connected to the viscoelastic 
foundation, which we briefly discuss in Subsection 11.11 We then briefly introduce the theory of 
Colombeau generalized functions which forms the framework for our work. Similar problems in- 
volving distributional and generalized solutions to Euler-BernouUi beam models have been studied 
in [4l [121 [ISl Uni mi- The development of the theory in the paper is divided into two parts. In 
Section [2] we consider the initial-boundary value problem ([3l)-(IC)-(BC) on the abstract level. We 
prove, in Theorem 12. 3[ an existence result for the abstract variational problem corresponding to 
(I3])-(IC)-(BC) and derive energy estimates which guarantee uniqueness and serve as a key 
tool in the analysis of Colombeau generalized solutions. In Section [31 we first show equivalence of 
the system (HI-© with the integro-differential equation ([3]), and apply the results from Section[5] 
to the original problem in establishing weak solutions, if the coefficients are in L°° . Afterwards we 
allow the coefficients to be more irregular, set up the problem and show existence and uniqueness 
of solutions in the space of generalized functions. 



1.1 The Euler-Bernoulli beam with viscoelastic foundation 

Consider the Euler-Bernoulli beam positioned on the viscoelastic foundation (cf. [5] for mechanical 
background) . One can write the differential equation of the transversal motion 

+ ^e[0,l],t>0, (5) 

where 

• A denotes the bending stiffness and is given by A{x) = EIi -(- H{x — xq)EI2- Here, the 
constant E is the modulus of elasticity, /i, /2, Ii ^ I2, are the moments of inertia that 
correspond to the two parts of the beam, and H is the Heaviside jump function; 

• R denotes the line density (i.e., mass per length) of the material and is of the form R{x) — 
Ro + H{x^xo)iRi-R2); 

• P{t) is the axial force, and is assumed to be of the form P{t) = Pq + PiS{t — ii), Po,Pi > 0; 

• g — gix,t) denotes the force terms coming from the foundation; 

• u = u{x,t) denotes the displacement; 

• h — h{x, t) is the prescribed external load (e.g. when describing moving load it is of the 
form h{x,t) — HqS{x — ct), Hq and c are constants). 

Since the beam is connected to the viscoelastic foundation there is a constitutive equation describ- 
ing relation between the force of foundation and the displacement of the beam. We use the Zener 
generalized model given by 

Dfuix, t) + u{x, t) = e Dfg{x, t) + g{x, t), (6) 

where < 6 < 1 and D" denotes the left Riemann-Liouville fractional derivative of order a with 
respect to <, defined by 

* T{l-a)dtJ^ (t-r)" 
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System (O-® is supplied with initial conditions 

u{x,0) = fi{x), dtu{x,Q) ^ f2{x), 

where /i and /2 are the initial displacement and the initial velocity. If /i (x) — /2 (x) — the only 
solution to (O-® should be u = g = 0. Also, the beam is considered to be fixed at both ends, 
hence boundary conditions take the form 

u{0, t) = t) = 0, 9^m(0, t) = d^u{l,t) = 0. 

By a change of variables 1 1-^ t via i(r) — ^ R{x)t the problem ([S])-® is transformed into the 
standard form given in ((I])-([2|). The function c in ([T]) equals A and therefore is of Heaviside type, 
and the function b is then given by b{x, t) = P{R{x)t) and its regularity properties depend on the 
assumptions on P and R. 

As we shall see in Section [31 standard functional analytic techniques reach as far as the follow- 
ing: boundedness of b together with sufficient (spatial Sobolev) regularity of the initial values /i, /2 
ensure existence of a unique solution u e i^(0, T; H§{{0, 1))) (in fact u S -L^(0, T; Hq{{0, 1)))) to 
dH]) with (IC) and (BC). However, the prominent case b = po + PiS{t — ti) is clearly not covered 
by such a result, so in order to allow for these stronger singularities one needs to go beyond 
distributional solutions. 

1.2 Basic spaces of generalized functions 

We shall set up and solve Equation ([3]) subject to the initial and boundary conditions (IC) and (BC) 
in an appropriate space of Colombeau generalized functions on the domain Xt ■= (0, 1) x (0, T) 
(with T > 0) as introduced in [3 and applied later on, e.g., also in [TTJ[T3]. As a few standard 
references for the general background concerning Colombeau algebras on arbitrary open subsets 
of R'' or on manifolds we may mention [S] IHl UHl [13 • 

We review the basic notions and facts about the kind of generalized functions we will employ 
below: we start with regularizing families {ue)e£{o,i] of smooth functions e H°°{Xt) (space of 
smooth functions on Xt all of whose derivatives belong to L^). We will often write {ue)e to mean 
(we)ee(o,i]- We consider the following subalgebras: 

Moderate families, denoted by £m.h°^(Xt)^ ^'"^ defined by the property 

VaeN^,3p>0: ||a"ue||L2(Xr) =0(£-''), as £ ^ 0. 

Null families, denoted by AfH°-={XT)j the families in £m,h°°{Xt) satisfying 

Wq>0:\\u,\\L2(^Xr)=Oie'') as e ^ 0. 

Thus moderateness requires estimates with at most polynomial divergence as £ — > 0, together 
with all derivatives, while null families vanish very rapidly as £ — ?> 0. We remark that null families 
in fact have all derivatives satisfy estimates of the same kind (cf. [H Proposition 3.4(ii)]). Thus 
null families form a differential ideal in the collection of moderate families and we may define the 
Colombeau algebra as the factor algebra 

Qh°^{Xt) = ^MM°°{Xt)/-^H°°(Xt)- 

A typical notation for the equivalence classes in Qh^(Xt) with representative {u^)^ will be 
[(ue)e]. Finally, the algebra Qh°°((o,i)) of generalized functions on the interval (0,1) is defined 
similarly and every element can be considered to be a member of Qh°°(Xt) well. 

We briefly recall a few technical remarks from [13l Subsection 1.2]: 

If {ue)e belongs to £m,H'^{Xt) have smoothness up to the boundary for every u^, i.e. 
Me S C°°([0, 1] X [0,T]) (which follows from Sobolev space properties on the Lipschitz domain 
Xt; cf. [1]) and therefore the restriction u\t=Q of a generalized function u S Gh°°(Xt) to t = is 
well-defined by Ue{-,G) € £m,H'^{{o,i))- 
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If u G Gh°°({q.i)) ^-iid in addition we have for some representative {ve)e of v that G Hq{{0, 1)), 
then ^^(O) = = and ^^^^^(O) = dxVe{l) — 0. In particular, 

v{0) = v{l) = and d^,v{0) = d^v{l) = 

holds in the sense of generalized numbers. 

Note that L^-estimates for parametrized families € H°°{Xt) always yield similar L°°- 
estimates concerning e-asymptotics (since H'^{Xt) C C°°(X^) C W'^-'°°{Xt)). 

The space H~°°{W^)^ i.e. distributions of finite order, is embedded (as a linear space) into 
QH°°(K.'i) by convolution regularization (cf. [3 ). This embedding renders H°°{U.'^) a subalgebra of 

Qh°° (R<i)- 

Certain generalized functions possess distribution aspects, namely we call u = [{u^)^] S Qh°° 
associated with the distribution w G V , notation w « w, if for some (hence any) representative 
(Mg)^ of u we have — > w in V, as e — > 0. 

2 Preparations: An abstract evolution problem in varia- 
tional form and the convolution- type operator L 

In this section we study an abstract background of equation ([3]) subject to the initial and boundary 
conditions (IC) and (BC) in terms of bilinear forms on arbitrary Hilbert spaces. First we shall 
repeat standard results and then extend them to a wider class of problems. We shall show existence 
of a unique solution, derive energy estimates, and analyze the particular form of the operator L 
appearing in ([3]). 

Let V and H be two separable Hilbert spaces, where V is densely embedded into H. Wc shall 
denote the norms in V and // by || • ||y and || ■ \\h respectively. If V denotes the dual of V, then 
V C H C V forms a Gelfand triple. In the sequel we shall also make use of the Hilbert spaces 
Ev := L'^{0,T;V) with the norm \\u\\e^ ■= {f^ \\u{t)\\l. dty^^ , and Eh := L^{0,T;H) with the 

norm ||u||_Ejf := {J^ \\u(t)\\'jj dt)^^'^ . Since, \\v\\h < C ■ \\v\\v, v eV, (without loss of generality we 
may assume that C = 1), it follows that HuUBfj < u G Ey, and Ey C Eh- The bilinear 

forms we shall deal with will be of the following type: 

Assumption 1. Let a(i, •,•), ao{t,-,-) and ai(i, •,•), t G [0,r], be (parametrized) families of 
continuous bilinear forms on V with 

a{t, u, v) — ao{t, u, v) + ai{t, u,v) Vu, w G V, 

such that the 'principal part' oq and the remainder ai satisfy the following conditions: 

(i) t I— ^ ao(i, u,w) is continuously differentiate [0,T] R, for all u,v G V; 

(a) flo is symmetric, i.e., ao(t, u,w) — ao{t,v,u), for all u,w G V; 

(Hi) there exist real constants A,/i > such that 

ao{t,u,u)> ^i\\u\\l,- X\\u\\%, VmG V^, VtG [0,r]; (7) 

(iv) t i-> ai{t,u,v) is continuous [0,T] — )• M, for all u,v G V; 

(v) there exists Ci > such that for all t G [0, T] and u,v gV , \a,i{t, u, v)\ < Ci||M||y HwHif ■ 

It follows from condition (i) that there exist nonnegative constants Co and Cg such that for all 
t G [0, T] and u, w G F, 

\aQ{t,u,v)\<Cr>\\u\\v\\v\\v and \a'^{t,u,v)\ < C'^\\u\\v\\v\\v . (8) 

where ao(i, u, v) := ^ao(i, u, v). 

It is shown in [3 Ch. XVIII, p. 558, Th. 1] (see also [W] Ch. HI, Sec. 8]) that the above 
conditions guarantee unique solvability of the abstract variational problem in the following sense: 
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Theorem 2.1. Let a{t, ■, •), t € [0,T], satisfy Assumption[l\ Let uq £ V, ui E H and f E Eh- 

Then there exists a unique u G Ey satisfying the regularity conditions 

u' = ^£Ev and u" ^ ^ G L\0,T;V') (9) 

(here time derivatives should be understood in distributional sense), and solving the abstract initial 
value problem 

{u"{t),v) +ait,u{t),v) ^ {f{t),v), yveV, for a.e. te{0,T) (10) 
u{0)=uo, m'(0)=ui. (11) 

(Note that (0) implies that u G C([0;T],1/) and u' G C([0, T], T^'). Hence it makes sense to 
evaluate u(0) G V and u'(0) G V' and claims that these equal uq and ui, respectively.) 

Remark 2.2. The precise meaning of ^TU\\ is the fohowing: V ip £ 2?((0,r)), 
((M"(i),u),<^)(x„^i,) + (a(t,'u(i),u),i^)(x)',-D) = ((/(i), u), 

or equivalently, 

rT j-T nT 

{u{t),v)f"{t)dt+ a{t,u{t),v)(p{t)dt = {f{t),v)ip{t) dt. 
Jo Jq 

The proof of this theorem proceeds by showing that u satisfies a priori (energy) estimates 
which immediately imply uniqueness of the solution, and then using the Galerkin approximation 
method to prove existence of a solution. An explicit form of the energy estimate for the abstract 
variational problem (|9l)- (|ll|) with precise dependence of all constants is derived in [13l Prop. 1.3] 
in the form 

lhWlly + h'WllH< (^Tho||^ + 11^^1111, + ^\/(r)||l,dr) -e*-^-, (12) 
^here Dt and Ft max{^^, ^i±^}. 



2.1 Existence of a solution to the abstract variational problem 

We shall now prove a similar result for a slightly modified abstract variational problem, which is 
to encompass our problem (|3]). Here in addition to the bilinear forms we shall consider "causal" 
operators L : L^{0,Ti;H) L^{0,Ti;H), VTi < T, which satisfy the estimate: 3 Cl > such 
that 

l|i"l|L2(0,Ti;H) < Cl\\u\\l'2(o.Ti-H), (13) 

where Cl is independent of Ti . 

Lemma 2.3. Let a{t, ■, ■), t G [0,r], satisfy Assumption[l[ Let fi £ V , f2 E H and h G Eh. Let 
L : Eh — > Eh satisfy Iil3\) . Then there exists a u £ Ey satisfying the regularity conditions 

, du ,, d^u 9 , 

u = — eEv and u" ^ — e 0, T; V') 
dt dt'' 

and solving the abstract initial value problem 

{u"{t),v) +a{t,u{t),v) + {Lu{t),v) ^ {h{t),v), \f v e V, for a.e. t £ {0,T), (14) 
uiO) = A, ^.'(0) = /2. (15) 

Moreover, we have u £ C{[0,T];V) and u' £ C{[0,T]; H). 
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Here we give a proof based on an iterative procedure and employing Theorem 12.11 and the 
energy estimate (|12p in each step. Notice that the precise meaning of (in distributional sense) 
is explained in Remark 12.21 

Proof. Let uq G Eh be arbitrarily chosen and consider the initial value problem for u in the sense 
of Remark [22] 

{u"{t),v)+a{t,u{t),v) + {Luo{t),v) ^ {h{t),v), Vu G for a.e. t G (0, T), (16) 
u(0)-/i, «'(0) = /2. 

By Theorem 12. II there exists a unique ui G Ey satisfying u[ G Ey, u" G L^(0,T; V), and solving 
P^ . Consider now ([T^ with Lui instead of Luq. As above, by Theorem 12. 1[ one obtains a 
unique solution U2 G Ey with G Ey and 1*2 G i^(0, T; F'). Repeating this procedure we obtain 
a sequence of functions {uk}k£N € E'y, satisfying uj. G ii'y, uj,' G L^(0,T; V), and solving the 
following problems: for each fc G N, 

«(i),w) +a(t,ufe(i),w) + (iufe_i(t),w) = (/i(t),v), yveV, for a.e. tG(0,T), 

Ufc(0)=/i, <(0)-/2. 

Also, for all fc G N, Ufc satisfies the energy estimate of type (fT2)) : 

IhfcWiPy + IKWIll^ < (^T II All'v + II/2III, + \\{h L«fe-i)(r)||l, dr) • e*'^-, 

where the constants Z?t and Ft are independent of k. We claim that {ufc}fcgN converges in Ey. 
To see this, we first note that ui — Uk solves 

{{ui - uk)"{t),v) + a{t,{ui - uk){t),v) + {L{ui-i - uk-i){t),v) =0, Vw G for a.e. te (0,r), 
{ui-uk){0) ^0, {ui-UkY {0)^0, 

and ui — Uk € Ey, with — m'^, G Ey and u" — w'^' G L^{0,T; V'). Moreover, the corresponding 
energy estimate is of the form 

Uui - UkKtWy + ||(u, - UkYml < e*-^^ • r ||i(«z-i - Uk-^)ir)rH dr. (17) 

Jo 

Thus, 

Wiui " ^^fc)(i)lly < e^-^" • / \\L{ui.i - Uk-i){T)\\l dr - e^'^- • ||L(wi_i - Ufc-i)|||„. 

Jo 

Integrating from to T and using assumption ()13p on L, one obtains 

||m; - Uk\\Ev < 7t||w;-1 - Uk-i\\EH, (18) 

where jt ClVTc 2 . Taking now Z = A; + 1 in ([18|) successively, yields 
||wfc+i - WfclUv < Trlki - ""olUff < 7tII"i " ""olUv; 

and hence 

k 

\\ui - UkWsv < 11"; - Ui-i\\ev + • • • + ||wfc+l - Uk\\Ev < X! '^tIWi - ^oWev 

i=l-l 

We may choose Ti < T such that < 1, hence X^i^oT'Ti converges. Note that t H> 7t is 
increasing. By abuse of notation we denote L'^{0,Ti;V) again by Ey. This further implies that 
{wfejfegN is a Cauchy sequence and hence convergent in Ey, say u := limfc^oo Uk- Similarly, one 
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can show convergence of u^, in Ey, i.e., existence of v := linifc_i.oo u'^, e Ey- In the distributional 
setting hmfe^.oo u'^, = u' , and therefore u' = v £ Ey (cf. 7, Ch. XVIII, p. 473, Prop. 6]). 
We also have to show that u solves equation p^ . Let ip e 2?((0,r)). Then 

{{h{t),v),p) = {{u'^{t),v),p) + {a{t,Uk{t),v),ip) + {{Luk-iit),v),p) 
= {{uk,v),ip") + {a{t,Uk{t),v),(p) + {{Luk-i{t),v),ip) 
^ ({u, v),ip") + {a{t, uit),v),p) + {{Luit), v),ip) 
= {{u'\t),v),ip) + {ait,u{t),v),ip) + {{Lu{t),v),ip). 

Here we used that (p" G I?((0,T)). Therefore u solves (ITi)) on the time interval [0,ri]. The initial 
conditions are satisfied by construction of u. 

It remains to extend this result on existence of a solution to the whole interval [0,T]. 

Since Ti is independent on the initial conditions, if T > Ti one needs at most ^ steps to reach 
convergence in Ey. In fact, one has to show regularity at the end point Ti of the interval [0,Ti] 
on which the solution exists, i.e., 

u{Ti) e V and u'(Ti) e H. 

To see this, it suffices to show that Uk ^ m in IV C{[0,Ti\;V) and u'j^ u' in Yh 
C{[0,Ti]; H). From ([T7)) and assumption ([T^ on L we obtain 

Uui - UkmWv < e^'-^^^Cl r \\{ui - Uk){T)\\ldT < e^-^-iCl ^ Um - Ufc)(r)|l^dT. 

Jo Jo 

Taking first the square root and then the supremum over all t £ [0,T] yields 

- Ufclln- < 7Ti||w; - Ukllvv 
Since < 1 this implies that {ufcjfceN is a Cauchy sequence in Yy. Similarly, 

\\{ui - u,ym% < e^-^-^d r \\{ui - Uk){T)\\l dr, 

Jo 

which upon taking the supremum gives 

\\{ui - UkYWrH < 7tiII"; - u/cIIfv: 

thus u'). — >■ u' in Yh (due to the already established convergence of Uk in Yy), and u'{Ti) € H. 
This proves the claim. □ 

2.2 Energy estimates 

In Section [3] we shall need a priori (energy) estimate for problem ([3]). In fact, for the verification 
of moderateness in the Colombeau setting it will be crucial to know all constants in the energy 
estimate precisely. Therefore, we shall now derive it. 

Proposition 2.4. Under the assumptions of Lemma \2.3l let u be a solution to the abstract vari- 
ational problem |J^[)-( f73)) . Then, for each t G [0,T], 

Mm'y + Wu'ml < (DTWhWl + - (ml + f WKrWudrX) ■ e*-^-, (19) 



V 

u • n 1 n Co + Ml+T) , Tj^ (CL+Ci+Cl Ci+2+A(1+T) i 

where v := mmjl, Dt '■= — ^^-^ — - and rr ■= max{— ^ — — -, — ■ — — 
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Proof. Setting v :— u'{t) in (jl4p we obtain (as an equality of integrable functions with respect to 
t) 

{u"{t),u'{t)) +a{t,u{t),u'{t)) + {Lu{t),u'{t)) = {h{t),u'{t)). 

Since a{t,u,v) = ao{t,u,v) + ai{t,u,v) and {u" (t) , u' {t)) = \f^{u'{t),u'{t)) = we 
have 

^Ik'WIIff = -2ao(<, ^i'W) - 2ai(<, u^, - 2{Lu{t), u'[t)) + 2{h{t), 
Integration from to ^i, for arbitrary < < T, gives 

IIMIh = -2/ aoit,u{t),u'it))dt-2 ai{t,u{t),u' {t)) dt 

{Lu{t),u'{t))dt + 2 [ {h{t),u'{t))dt. 
Jo 

Note that ^aQ{t,u{t),u{t)) — a'Q{t,u{t),u{t)) + aQ{t,u'{t),u{t)) + aQ{t,u{t),u'{t)) and hence, by 
Assumption [T] (ii) , 2aQ{t,u'{t),u{t)) = -^aQ{t,u{t),u{t)) — a'Q{t,u{t),u{t)). This yields 



LHS := \\u'{h)\\l + aoih,u{ti),u{ti)) = ||/2|||j + ao(0, u(0), m(0)) - / a'o{t,u{t),u{t)) dt 

Jo 

ai{t,u{t),u'{t))dt^2 I {Lu[t),u'{t))dt + 2 [ {h{t) , u' [t)) dt RH S . (20) 



Further, by ([5]), Assumption [T] (v), the Cauchy-Schwartz inequality, the inequality 2ah < + 5^, 
and the assumption (|13p on L we have 



\RHs\<\\hrH+c4umi+c'^ / \\u{t)\\idt 

Jo 

+ 2Ci r \\u{t)\\v\\u'{t)\\Hdt + 2 r \\Lu(t)\\H\\u'{t)\\Hdt + 2 I \\h{t)\\H\\u'{t)\\H dt 



< IIMIlf + Co||/i|py + (C^ + Ci) / \\u{t)\\l.dt 

Jo 

+ (Ci+2) f'\\umldt + \\Lu\\l.(,^,^.^H)+ r \\h{t)\\ldt 
JO Jo 

< \\f2fH + Co\\fifv + {C^ + Ci+CL) r\\u{t)\\ldt 

Jo 

+ (Ci+2) r \\u'{t)\\]jdt+ r \\h{t)\\]jdt. 



Further, it follows from ([7]) that 

LHS = \\u{h)\\% + ao(ii, w(ii), w(ii)) > h'(ii)lll/ + ^J'\Htl)\\l - \\Hti)\\]j, 
and therefore ((20| yields 







+ (C^ + Ci+Cl) /*'||u(t)||2,dt + (Ci+2) / h'(t)|l^di. 



As shown in [13] we have that \\u{t)\\]j < (1 +<)(||/i||y + ||u'(s)|||^ ds), hence 



where v := min{l,^}, Dt := and Ft := max{ , ^1+2+^^(1+^) }. The clahn 

now fohows from Gronwall's lemma. □ 

As a consequence of Proposition l2.4[ one also has miiqueness of the solution in Lemma [2.31 

Theorem 2.5. Under the assumptions of Lemma \2.3\ there exists a unique u G Ey satisfying the 
regularity conditions u' € Ey and u" e L^{0,T; V'), and solving the abstract initial value problem 
(T^-dMl- Moreover, ueC{[0,T];V) and u' e C{[0,T]; H). 

Proof. Since existence of a solution is proved in Lemma 12. 3i it remains to show uniqueness part 
of the theorem. Thus, let u and w be solutions to the abstract initial value problem P^ - ([T5|) . 
satisfying the regularity conditions u' , w' G Ey and u", w" G L^{0, T; V'). Then u — w is a solution 
to the homogeneous abstract problem with vanishing initial data 

{{u-w"{t),v) +a(t, {u-w){t),v) + {L{u~ w){t),v) ^ 0, Vw G V, for a.e. t e (0,r), 
(m-u;)(0)=0, (w - u))'(0) = 0. 

Moreover, according to Proposition 12. 4[ u — w satisfies the energy estimates ([TO)) with /i = /2 = 
h = 0. This implies uniqueness of the solution. □ 

2.3 Basic properties of the operator L 

In this subsection we analyze our particular form of the operator L, relevant to the problem 
described in the Introduction. Therefore, we consider an operator of convolution type and seek 
for conditions which guarantee estimate p^ . 

Lemma 2.6. Let I e Lf^^{R) with supp/ C [0, 00). Then for all Ti e [0,T], the operator L 
defined by Lu{x,t) :— l{s)u{x,t — s) ds maps L^{0,Ti] H) into itself, and US]) holds with Cl = 

\\l\\L^Q,T) ■ T. 

Remark 2.7. We may think of u being extended by outside [0,T] to a function in L^(IR;_ff), 
and then identify Lu with / *t u. 

Proof. Integration of < /J - s)\\\u{s)\\Hds from to Tj, < Ti < T, yields 




< IKI|l2(o,T) • T ■ ||u||l2(o^Ti://), 

where we have used the support property of I, Minkowski's inequality for integrals (c.f |^ p. 194]), 
and the Cauchy- Schwartz inequality. □ 

In the following lemma we discuss a regularization of L, which will be used in Section [X^ 

Lemma 2.8. Let I e L\^^{R) with suppZ C [0,oo). Let p £ 2?(M) be a mollifier fsuppp C -Bi(O), 
J p = 1). Define Pe{t) := ^ePilet), with 7^ > and ^ 00 as e ^ 0, l^ :— I * p^ and 
L^u{t) {l^ *t u){t), for u G Er. Then Vp G [l,oo), l^ G ^LW (^"^^ l^ I in i^W- 
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Proof. Let K he a, compact subset of M. Then 

i/p 



\k\\LP{K) = \\l *t PeWLPiK) = (^J \j 



l{T)p^{t~T) dT\Pdt 



p s 1/p / „ / „ \ P \ 1/P 



< / / \Kr)\\Mt~r)\dr dt < / / |Ur)||pe(t - r)| dr dt \ 

\Jk \J-oo / J \Jk \Jk+Bi{0) ) ) 

<f {[ \l{T)npS-rrdt\'' dT= [ \l{T)\([\pS-TWdt\'\T 

JK+Bi(Q) \Jk J Jk+Bi(0) \Jk J 

l^('^)ll|Pe||LJ'(Bi(0)) dT = ll^llLi(^+-Bi(0))ll/'£lli''(Si(0)) 

K+-Bi(0) 

1-i 

= IKIlLi(if+Bi(0)) • 7e " ■ IIpIIlp(Bi(0))- 

where the second inequahty fohows from the support properties of I and p {t — t ^ Bi{0),t G K 
imphes r G K + Bi(0)), while for the third inequahty we used Minkowski's inequahty for integrals. 
Further, we shaU show that ^ / in Ll^^{M.). Let K CC R. We claim that \le — Z| — > 0, as 
e -> 0. Indeed, 

I / l{t~ s)p,{s)ds- l{t) ■ / peis)ds\dt^ / I / {l{t~ s)~l{t))peis)ds\dt 



{l{t )-l{t))p{T)dT\dt < / / \l{t )~l{t)\\p{T)\dTdt 

k'Jr le JkJr 7e 

= / \pir)\ I \l{t--)-l{t)\dtdT. 

By [8l Prop. 8.5], we have that ||/(- — —) — ^IIli(a') — >■ 0, as e — > and therefore the integrand 
converges to pointwise almost everywhere. Since it is also bounded by 2|p(T)|||Z||ii(A) € 
Lebesgue's dominated convergence theorem implies the result. □ 



3 Weak and generalized solutions of the model equations 

We now come back to the problem |T])-([2])-(IC)-(BC) or ([3])-(IC)-(BC), and hence need to provide 
assumptions which guarantee that it can be interpreted in the form (|14p . in order to the apply 
results obtained above. For that purpose we need to prescribe the regularity of the functions c 
and h which appear in Q. In Section [5^ we shall use these results on the level of representatives 
to prove existence of solutions in the Colombeau generalized setting. 

Thus, let H :— i^(0, 1) with the standard scalar product = j^u{x)v{x)dx and i^- 

norm denoted by || • \\h- Let V be the Sobolev space i?Q((0,l)), which is the completion of 
the space of compactly supported smooth functions C^((0, 1)) with respect to the norm ||u||2 = 
(ELo (and inner product {u,v) ^ Y1=o{'^^^\v^^^))- Then V = i?-2((0,l)), which 

consists of distributional derivatives up to second order of functions in L^(0, 1), and V ^ H ^ 
V' forms a Gelfand triple. With this choice of spaces H and V we also have that Ey = 
L2(0,T;7J2((0,1))) and Eh = L^ {{0,1) x (0,T)). 

Let 

ceL°" (0,1) and real, & e C([0, T]; L°°(0, 1)), (21) 
and suppose that there exist constants ci > cq > such that 

< Co < c{x) < ci, for almost every x. (22) 

For t £ [0,r] we define the bilinear forms a{t, •, •), ao(t, •, •) and ai{t, ■,-)onVxV by 

ao{t,u,v) — {c{x) d'^u,d'^v), ai{t,u,v) — {b{x,t) d^u,v) , u,v&V, (23) 
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and 

a{t, u, v) = ao(t, u, v) + ai{t, u, v). (24) 

Properties (Hi}, ([2^ imply that qq, ai defined as in (|23p satisfy tlie conditions of Assumption [T] 
(cf. [13, proof of Th. 2.2]). The specific form of the operator L is designed to achieve equivalence 
of the system ©-(HI) with the equation which we show in the sequel. 

Let denote the space of Schwartz' distributions supported in [0, oo). It is known (cf. \18\ ) 
that for given z S S'^ there is a unique y e S'^ such that -Df z + z — ODfy + y. Moreover, it is 
given hy y = Lz, where L is linear convolution operator acting on S'j^ as 

f^^^ ^^S' (25) 



The following lemma extends the operator L to the space Eh- 

Lemma 3.1. Let L : S'_^_ — > S'_^_ he defined as in 1125]}. Then L induces a continuous operator 
L = Id+La on Eh, where La corresponds to convolution in time variable with a function G 

Proof. Recall that for the Mittag-Leffler function ea{t. A), defined by 

e it X)-Y^~^ 

^"(*'^)-l.r(.. + i)' 

we have that £{ea{t, X)){s) = f^, € C°°((0,oo)) nC([0,oo)) and e'„ G C°°((0,oo)) n 
iL([0,oo)) (cf. m)- Also, 

( it) ^C-'(l + ^^-^l it) ^ S{t) +(l-l]e'Jt, I) 5{t) + L{t). 



Let ue Eh. Then 

(^^) (•) *t "(^, •) = "(^, •) + - l) e'a * nix, .) (26) 

is an element in L^{0,T) for almost all x (use Fubini's theorem, e'^ £ L^{0,T) and L^ * L^ C L^ 
(cf. [H]))- Extend this to a measurable function on (0,1) x (0,r), denoted by Lu. By Young's 
inequality we have 

\\{Lu){x,-)\\l^o,t) < \\u{x,-)\\l^o,t) + 1^- M\\e'c*u{x,-)\\L^o,T) 

< Mx, •)IIl2(o.t) + - l|l|eallLi(o.T)l|u(a;, ■)\\l^o.t), 

hence, 

\\Lu\\e, < (1 + - m<MLHo.T))\\u\\E,. (27) 
Thus, Lu G Eh and L is continuous on Eh. □ 
We may write 

Lu :— (Id +La)u — I *t u — {S + la) *t u with LaU :— la *t u, la ■= {-p: — l)eQ(i, -), (28) 



and therefore the model system ([I])-© is equivalent to Equation ([3]). 
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3.1 Weak solutions for L°° coefficients 

Now we are in a position to apply the abstract results from the previous section to the original 
problem. 

Theorem 3.2. Let b and c be as in i21]) and \2'2\) . Let the bilinear form a{t, ■, ■), t £ [0,T], be 
defined by H^j and (2^, and the operator L as in (23\). Let fi G Hq{{0,1)), /2 € L^{0,1) and 
h e L^((0, 1) X (0,T)). Then there exists a unique u £ L^{0,T; Hq{0, 1)) satisfying 

u' = ^e L2(0,T;i/2(0,l)), = ^ e L'{0,T;H~\0,1)), (29) 

and solving the initial value problem 

{u"{t),v)+a{t,u{t),v) + {Lu{t),v)^0, Vu G iJo^((0, 1)), i G (0, T), (30) 
u(0)=/i, u'{0) = f2. (31) 

(Note that, as in the abstract version, since \29\} implies u £ C([0, T], iJp ((0, 1))) and u' € 
C{[0,T],H-'^{{0,1))) it makes sense to evaluate u{0) G iJ^((0, 1)) and li'(O) G ff"^((0,l)) and 
i31\) claims that these equal fi and f^, respectively.) 

Proof. We may apply Lemma [^751 because the bilinear form a and the operator L satisfy Assump- 
tion[T]and condition p3)) . The latter is true according to (|27)) with Cl = (1 + 1^ — l|||ea||Li(o,r)) = 
1 + II'q||l2(o,t)- As noted earlier, the bilinear forms a, oq and a\ are as in [131 (20) and (21)]. 
Moreover, it follows as in the proof of [13j Theorem 2.2] that a satisfies Assumption [1] with 

Co ||c||l=°(04), Cq := 0, Ci ||6||loo((o^i)x(o,t)), M A := C1/2 • co, (32) 

where C1/2 is corresoponding constant form Ehrling's lemma. □ 

We briefiy recall two facts about the solution u obtained in Theorem l3.2l (as noted similarly in 
[El Section 2]): 

(i) Since u{.,t) e i/Q((0, 1)) for ah t e [0,T] and i/o((0, 1)) is continuously embedded in 
{v e Ci([0, 1]) : i;(0,t) = w(l,i) = 0,a^w(0,t) = d^v{\,t) = 0} ([19, Corollary 6.2]) the solution u 
automatically satisfies the boundary conditions. 

(ii) The properties in (gH) imply that u belongs to Ci([0, T], i7-2((0, 1))) n L2((o, T) x (0,1)), 
which is a subspace of I?'((0, 1) x (0, T)). Thus in case of smooth coefficients b and c we obtain a 
distributional solution to the "integro-differential" equation 

diu + dl{cdl)u + bd'^u + Utu^ h. 

3.2 Colombeau generafized solutions 

We will prove unique solvability of Equation ^ (or equivalently, of Equations ([I])-©) with (IC) 
and (BC) for u G Gh°°{Xt) when b,c, fi, f2,g and h are Colombeau generalized functions, where 
Xt:= (0,1)x (0,r). 

In more detail, we find a unique solution u G Qh°°(Xt) the equation 

dfu + Q{t,x,dx)u + Lu = h, on Xt 

with initial conditions 

u\t=a = /i G Gh^({o,i)), dtu\t=o = /2 e Gh°°({o,i)) 
and boundary conditions 
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Here Q is a partial differential operator on Gh=°{Xt) with generalized functions as coefficients, 
defined by its action on representatives in the form 

Furthermore, the operator L corresponds to convolution on the level of representatives with reg- 
ularizations of / as given in Lemma [27 



{Ue)e ^ {Is *t Ue{t))^ -. (L^Ue 

where 1^ — I * Pe, with introduced in Lemma [ 



Lemma 3.3. (i) If I e £L(M) and I isC^ in {0,oo) then L is a continuous operator on H°° {Xt)- 
Thus ^ {Lu^)^ defines a linear map on Gh=°{Xt)- 

(a) If I G Lj^^(M.) then Ve € (Oj 1] the operator is continuous on II°°(Xt) and (u^)^ i— > 
{Lue)e defines a linear map on Qh=°(Xt)- 

Proof, (i) From Lemma 12.61 with H = i^(0, 1) we have that L is continuous on L'^{Xt) with 
operator norm ||i||op < T ■ \\1\\l^(q^t)- 

Let u G H°°{Xt) and Lu{x, t) = /J l{s)u{x, t — s) ds. We have to show that all derivatives of 
Lu with respect to both x and t are in IF'(Xt)- 

• di.Lu{x,t) = J*l{s)di.u{x,t-s)ds, and hence, \\dlLu\\L2(^XT) < \\l\\L^o,T)\\diu\\L2(^XT)- 

• d^dl.Lu = d^L{dl.u), and since the estimates for L{dl.u) are known it suffices to consider 
only terms d^Lu. For the first order derivative we have 

dtLu{x,t) — l{t)u{x,0) + / l{s)dtu{x,t — s) ds 

Jo 

and therefore 

WdtLuU^^Xr) < ll^llL^(o,T)ll"(-,0)|U2(oa) + T. |l/|U.(„,T,)|15tw|U2(;,,) 
< \\1\\lHo,t){\\u\\h'^{Xt) + T ■ \\u\\h^(Xt)), 

where we have used the fact that Tr : II°°{Xt) — t- -ff°°((0, 1)), u i— >■ w(-,0) is continuous, 
and more precisely, Tr : H'^iXr) ^ iJ^^^lO, 1)) with estimates \\dld^u{-,0)\\L2^o,i) < 
l|w||ffm(Xr): = m{k,l). 

Higher order derivatives involve terms l^^\t)dlu{x, 0), . . . , /g l{s)d^u{x, t — s) ds, which can 
be estimated as above. 



(ii) From Lemma[278]it follows that l^ e -^focW: and ||/£|1l2(o,t) < li ■ I1^I1li(o,t)I1/o|U2(o,t)- 
From Lemma 12.61 we know that is continuous Xt — >■ Xt, with ||ie||op < T ■ ||^e||L2(o,T) ^ 

T -ji • ||^|1li(o.t)||pI|l2(o.t)i which is moderate. We can now proceed as in (i) to produce estimates 

of ||ieWe||i/^(XT): e always replacing \\1\\l^o,t) by ji ■ ||'||li(o,t)IIpIIl2(o,t) factors. Since 

'■'H°°{Xt)- 



7£ < e ^ it follows that {L^u^)^ £ Sff^rxT)- ^ 



Remark 3.4. The function I as defined in ^ belongs to Lf^^{R), if a > 1/2, and to Lj^^iR), if 
< 1/2 (which follows from the explicit form of e^(t, |)). This means that in case a > 1/2 we 
could in fact define the operator L without regularization of I. 

As in the classical case we also have to impose a condition to ensure compatibility of initial 
with boundary values, namely (as equation in generalized numbers) 

/i(0) = /i(l) = 0. (33) 
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Note that if /i G Gh=°({o.i)) satisfies then there is some representative {fi.e)s of /i with the 
property /i.^ G Hq{{0, 1)) for aU e G (0, 1) (cf. the discussion right below Equation (28) in [T3]). 

Motivated by condition ([2^ above on the bending stiffness we assume the following about c: 
There exist real constants ci > cq > such that c e Gh°^ [o.i)) possesses a representative {ce)e 
satisfying 

< Cq < c,{x) < ci Vx e (0,l),Ve e (0,1]. (34) 

(Hence any other representative of c has upper and lower bounds of the same type.) 

As in many evolution-type problems with Colombeau generalized functions we also need the 
standard assumption that b is of L°°-log-type (similar to [H]), which means that for some (hence 
any) representative (5e)e of b there exist N Gf^ and £o G (0, 1] such that 

||&e||L~(x^) <A^-log(J), 0<e<eo. (35) 

It has been noted already in fW, Proposition 1.5] that log- type regularizations of distributions are 
obtained in a straight-forward way by convolution with logarithmically scaled moUifiers. 

Theorem 3.5. Let b G Gh°^(Xt) of L°° -log-type and c G Gh°°{o.i)) satisfy (3^. Let 7^ = 
O(logi). For any fi G Gh-^hoa)) satisfying f2 G Gh-^((o,i)), h G Gh°°{Xt) and I G 

Gh°°((o,1))7 there is a unique solution u G Gh°°(Xt) initial-boundary value problem 

dfU -\- Q{t, X, dx)u + Lu = h, 
u\t=o = /i, dtu\t=o = /2, 

wU=0 = wU=l = 0, dxU\x=0 = da;U\x=l = 0. 

Proof. Thanks to the preparations a considerable part of the proof may be adapted from the 
corresponding proof in [13, Theorem 3.1]. Therefore we give details only for the first part and 
sketch the procedure from there on. 

Existence: We choose representatives (&e)e of b and {c^)e of c satisfying (p4l) and (|35|) . 
Further let (/ie)e, (/2e)e, ('e)ej and (/ie)e be representatives of /i,/2, I, and g, respectively, where 
we may assume /i,^ G H^HO, 1)) for all e G (0, 1) (cf. (|33l) ). 

For every e G (0, 1] Theorem l3 . 2l Drovides us with a unique solution G H^{{0, T), i?o((0, l)))n 

ij2((o,r),i?-2((o,i)))to 

PgWe :— dfUe -\- Qe{ii dx)Ue + L^Ue — On Xt, (36) 

'"e|t=0 = /l£, dtUe\t=(i ^ f2e- 

In particular, we have G CH[0, T], i?-2((-o, 1))) n C([0, T], iJ^((0, 1))). 
Proposition [23] implies the energy estimate 



\We{tWm + \\<mh < {D^T Wfiefm + \\f2dh + f \\Ut)\\1. dr) ■ exp(tFI), (37) 

where with some N we have as e — > 
D't = (llcelk- + A(l + T))/mm{n, 1) = 0{\\c,\\l^) - 0(1) (38) 

= -^--{\M^-+C^M^ +2 + Hl+T)} ^ ^ ii^^i,^^^ ^ 0(log(s--)), (39) 

mm(/i, 1) 



since /i and A are independent of e, and Cl.e = 0(log -) (cf. ([32])). 

By moderateness of the initial data /i^, f2e and of the right-hand side the inequality (|37p 
thus implies that there exists M such that for small e > we have 

+ lia.T/elli.fX.) + l|5^^e|li.(X.) + \\dtU,\\h(^^) = 0(£~'^), s ^ 0. (40) 

From here on the remaining chain of arguments proceeds along the lines of the proof in [131 
Theorem 3.1]. We only indicate a few key points requiring certain adaptions. 
The goal is to prove the following properties: 
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1. ) For every e € (0, 1] we have £ H^{Xt) C C°°(Xt). 

2. ) Moderateness, i.e. for all ^, fc G N there is some Af G N such that for small £ > 

Note that (gO]) already yields for {I, k) G {(0, 0), (1, 0), (0, 1), (0, 2)}. 



Proof of 1.) Differentiating (|36l) (considered as an equation in I?'((0, 1) x (0, T))) with respect to 
t we obtain 

Pe{dtUe) = dth^ - dtbe{x,t)dlu^ - le{t)fl,e =■ he, 

where we used dt{LeU,) = L^idtUe) + le{t)ue{0). We have G H^{{0,T), L'^{0, 1)) since dtK G 
H^{Xt), k e i7°°((0,T)), /i.e G i/°°((0,l)), at6e(a;,0 e H"°{Xt) C M^°°^°°(Xt) and ^^^^ g 
H^{{0,T),L^{0, 1)). Furthermore, since depends smoothly on i as a differential operator in x 
and Ue(0) = fi,e e i?°°((0, 1)) we have 

(5tU,)(-,0)-/2,e-:/l.eei^°°((0,l)), 
(5t 0) = K{; 0) - ge(We(-, 0)) - L,U,{-,0) = 0) - (Qe + /2,e G i?°"((0, 1)). 

Hence dtu^ satisfies an initial value problem for the partial differential operator as in p6l) 
with initial data /i^e, f2,e and right-hand side instead. However, this time we have to use 
V = i/^((0,l)) (replacing H^{{0,1))) and H = L^(0,1) in the abstract setting, which still can 
serve to define a Gelfand triple V ^ H ^ V' (cf. [El Theorem 17.4(b)]) and thus allows for 
application of Lemma [2731 and the energy estimate ([T9| (with precisely the same constants). 

Therefore we obtain dtUe G H^{[0,T], H'^{{0, 1))), i.e. G H'^{{0,T), H'^{{0, 1))) and from the 
variants of (|57)) (with exactly the same constants and and (pU| with 9tWe in place of 
that for some M we have 

+ + + = 0(e-*^) (e ^ 0). (41) 



Thus we have proved {Ti^k I with (/, k) = (2, 0), (1, 1), (1, 2) in addition to those obtained from (j40p 
directly. 

The remaining part of the proof of property 1.) requires the exact same kind of adaptions in 
the corresponding parts in Step 1 of the proof of 13, Th. 3.1] and we skip its details here. In 
particular, along the way one also obtains that 

(Ti^k) holds for derivatives of arbitrary I and k < 2. 

Proof of 2.) From the estimates achieved in proving 1.) and equation (j36]) we deduce that 

fee := a|(Ce d^Ue) = - dlue - dfue - L^Ue 

satisfies for alH G N with some Ni an estimate 

||5^fce|U.(x.) = O(£-^0 (£^0). (42) 

Here we are again in the same situation as in Step 2 of the proof of [13l Theorem 3.1], where 
now ke plays the role of /i^ there. Skipping again details of completely analogous arguments we 
arrive at the conclusion that the class of {ue)e defines a solution to the initial value problem. 

Moreover, we have by construction that Ui;{t) G Hq{{0,1)) for all t G [0, T], hence u{0,t) — 
u{l,t) — and dxu{0,t) — dxu{l,t) = and thus u also satisfies the boundary conditions. 

Uniqueness: If u = [(«£)£] satisfies initial-boundary value problem with zero initial values 
and right-hand side, then we have for all g > 

ll/i^ell =0(e«), ||/2.e|| =0(6''), \\K\\L2^Xr) ^ 0{e'^) as £ ^ 0. 

Therefore the energy estimate p7l) in combination with (|38p - ([5^ imply for all g > an estimate 

\K\\mxr)^0{e^) (e^O), 

from which we conclude that (Me)^ G J^h=°(Xt)^ ^-^-j = 0. □ 
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